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Abstract 

It is shown that the Euler system of hydrodynamic equations for inviscid 
barotropic fluid for density and velocity is not a complete system of dynamic 
equations for the inviscicd barotropic fluid. It is only a closed subsystem of 
four dynamic equation. The complete system of dynamic equation consists of 
seven dynamic equations for seven dependent variables: density, velocity and 
labeling (Lagrangian coordinates, considered as dependent variables). Solu- 
tion of the Cauchy problem for the Euler subsystem is unique. Solution of 
the Cauchy problem for the complete hydrodynamic system, containing seven 
equations, is unique only for irrotational flows. For vortical flows solution of 
the Cauchy problem is not unique. The reason of the nonuniqueness is an 
interfusion, which cannot be taken into account properly in the framework of 
hydrodynamics. There are some arguments in favour of connection between 
interfusion and turbulence. 

1 Introduction 

Describing a fluid flow, one labels particles of the fluid by means of Lagrangian coor- 
dinates £ = {61,62; £3}- O ne supposes, that any Lagrangian coordinates £ label the 
same fluid particle all the time. Applying laws of Newtonian dynamics to any fluid 
particle, one obtains hydrodynamic equations for the fluid flow in the Lagrangian 
representation (in the Lagrangian coordinates). The Lagrangian representation is 
sensitive to the correct labeling of the fluid particles in the sense, that there are situ- 
ations, when in different time moments the same Lagrangian coordinates £ describe 
different fluid particles. For instance, if two like gas beams, consisting of nonin- 
teracting molecules, pass one through another, the particle labeling changes after 
"collision" of the two beams. The picture is shown in the figure, which describes 
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world lines of gas particles in the space-time. The solid lines shows gas particle with 
the same labeling whereas dashed lines show world lines of real gas molecules. 

In this example a violation of the fluid particle labeling after "collision" is evident. 
The stream lines, represented by solid lines, do not describe motion of real fluid 
particles. 

Is it important? Can we observe stream lines of a fluid? Let us imagine that 
we introduce several flecks of dust in one of gas flows and follow their motion. We 
suppose that the size of flecks is larger, than the size of gas molecules, and any 
fleck interacts with many gas molecules. Then we may think, that any fleck of dust 
moves along stream line of the fluid. As far as collision of any fleck of dust with gas 
molecules is random, the flecks of dust after "collision" appear in both gas beams, 
although before the "collision" they were placed only in one of them. On the other 
hand, it is generally assumed, that flecks of dust moves together with the gas, and 
any fleck moves along the stream line of the fluid. Observation of flecks of powder 
is a usual method of the stream lines investigation. 

The hydrodynamic description of a fluid is valid at the supposition, that one 
stream of a fluid can penetrate into the other one only to the depth of the mean 
length of collision path. In this case the interfusion of different streams of a fluid 
will be infinitesimal. However, such a small interfusion will take place, and this 
interfusion may appear to be essential for the shape of stream lines. 

Such a physical phenomenon as turbulence can be discovered only, if one traces 
the irregular behavior of stream lines. In other words, for observation of turbulence 
a displacement of fluid particles is important, but not only their velocities. The 
velocities are important only as a source of displacement. 

A motion of the inviscid barotropic fluid is described by the Euler equations 



_ + (vV)v = —Vp, P = p(p) = P 2 -t^ (1.1) 

^ + V(pv) = (1.2) 

where p = p(t, x) is the fluid density, v = v (i, x) is the fluid velocity, p — p(p) = 
p 2 dE (p) I dp is the pressure and E (p) is the fluid internal energy per unit mass. 

Stream lines are described by a system of ordinary differential equations. One 
supposes that fluid particles move along the stream lines x = x (t), which are defined 
by the equation 

^ = v(f,x) (1.3) 

where v = v (t, x) is a solution of the Euler equations (11.11) . (11.21) . 

The Euler system (11. ip . (I1.2p is a closed system of differential equations, which 
may be solved independently of equations (11.31) . The Euler system is a system of 
nonlinear partial differential equations. It is difficult for solution. The system of 
ordinary differential equations (II. 3p is simpler, than the Euler system. Besides, it 
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can be solved only after solution of the Euler system (11. ip . (11.21) . It is a reason, 
why researchers investigate mainly the Euler system. The system of equations (II. 3p 
for the stream lines remains usually outside the region of consideration as some 
triviality. 

However, the Euler system (11.11) . ( 11. 2ft and equations (11.31) are dynamic equations 
of one dynamic system, and they should be considered together. This dynamic 
system will be referred as the complete hydrodynamic system. The Euler equations 
(11.11) . (11.21) are not dynamic equations of a wholesome dynamic system, because the 
cannot be deduced from a variational principle, whereas dynamic equations ( 11.11) . 
( 11.21) . ( 11. 3p can. It seems, that dynamic equations do not influence on the solution 
of dynamic equations (11.11) . (11. 2p . In reality, it is true only for irrotational flows. In 
vortical flows a situation changes. In the vortical flows an interfusion appears. The 
interfusion is not so large as in figure with colliding gas beams. This interfusion is 
infinitesimal. It is conditioned by different velocities of adjacent fluid volumes. This 
infinitesimal interfusion influences on the shape of stream lines and on the labeling 
of the fluid particles, although it does not influence on quantities p, v, which are 
solution of the Euler system (jl.ip . (11.21) . 

In the case of the Euler representation of hydrodynamic equations, the fluid par- 
ticle labeling is not necessary. The four hydrodynamic Euler equations are obtained 
as a result of the conservation laws of the energy and momentum. However, in 
this case the conservation law of the angular momentum is not used, and one can- 
not be sure, that the system of four Euler equations for barotropic fluid describes 
completely rotational degrees of freedom of molecules and those of fluid particles. 
The rotational degrees of freedom may be essential in turbulent flows, where close 
consideration of rotational degrees of freedom may appear to be essential. 

In this paper we try to take into account interfusion, which appears in the rota- 
tional flows of the barotropic fluid. Influence of interfusion manifests itself in labeling 
of the fluid particles by means of the Lagrangian coordinates £ = {£1,^25 £3}- This 
change of labeling is not important for solution of Euler equations (jl.ip . (11.21) . which 
form a closed system of differential equations. However, this change of labeling may 
appear to be important for such physical phenomena, where the shape of stream 
lines is essential (such as turbulence). 

We consider connection between the labeling and the interfusion on the for- 
mal mathematical level. We shall consider the Euler dynamic equation for inviscid 
barotropic fluid 

We are interesting in the question, whether the labeling, generated by the equa- 
tion (ll.3p is an unique possible way of labeling. To solve this problem, we shall 
consider the Lagrangian coordinates £ to be dependent dynamic variables. The 
Eulerian coordinates t, x are considered to be independent dynamic variables. 

Thus, we consider dynamic system, described by seven dependent dynamic vari- 
ables p, v, which are functions of four independent variables t, x. Note that system 
of dynamic equations (11.11) . (I1.2p is a closed system of dynamic equations. However, 
the dynamic system, described, by four dependent variables p, v is not a wholesome 
dynamic system in the sense, that dynamic equations (11.11) . (11.21) cannot be obtained 
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from some a variational principle. 

To obtain dynamic equations (11. ip . (11. 2p from the variational principle |5], one 
needs to add so-called Lin constraints [3]. This conditions have the form 

<9o£+(W)£ = (1.4) 

It is easy to see, that characteristics of the linear differential equation (11.41) 

dit doc cIjOlj d/0c , >, 

T = ~T~ = — = ~Y L5 

1 V 1 V z V 6 

coincide with the equation (jl.3p . Vice versa, any integral of the equation system (II. 3p 
is a solution of the equation (II .4p . The Lin constraints (11.4j) are interesting in the 
relation, that independent dynamic variables in (11.41) are the same, as in the dynamic 
equations (II. ip . (I1.2p . Hence, dynamic equations (11.11) . (I1.2p and dynamic equations 
(11.41) may be considered as dynamic equations of one dynamic system. It is rather 
difficult to consider system of equations (II. ip . (11.21) . (11.31) as a dynamic equations 
of a dynamic system, because independent variables are different in equations (11.11) . 
(1L2D and (Q. 



2 Generalized stream function 

Let us note that the quantities £ may be considered to be the generalized stream 
function (GSF), because £ have two main properties of the stream function. 

1. GSF £ labels stream lines of a fluid. 

2. Some combinations of the first derivatives of any £ satisfy the continuity 
equation identically. 

d k j k = 0, J k = ^, d k = J^, k = 0,1,2,3 (2.1) 

where j h = {j°, j l ,j 2 , j 3 } = {p, pv} is the 4-vector of flux. Here and in what follows, 
a summation over two repeated indices is produced (0-3) for Latin indices and (1-3) 
for Greek ones. The Jacobian determinant J = J^/ x 

J$/x = J fc, fc ) = g ^ = det I I ' = l,k = 0,1,2,3 

(2-2) 

is considered to be a four-linear function of £ lk . The quantity £ is the temporal 
Lagrangian coordinate, which appears to be fictitious in expressions for the flux 
4-vector j k 

_ 9 (f 1, $2, 6$) 1 _ ,1 _ 9 (f 1, f 2 , £3) 



J = Q /_1 _9 _ 3 N » P» =J 



d (x l , x 2 , x^Y ^(t,x 2 ,x 3 ) , 

2 -2 _ d(£i,£ 2 ,£ 3 ) m .3 _ -3 _ 9 €2, €3) (9 o\ 

PV ~ 3 ~ d(t,x\x*)> PV ~ 3 ~ d(t,x\x 2 ) [2 ' 6) 
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A use of Jacobians in the description of the ideal fluid goes up to Clebsch [H 
[2j Clebsch, 1857,1859], who used Jacobians in the expanded form. It was rather 
bulky. We use a more rational designations, when the 4-fl.ux and other essential 
dynamic quantities are presented in the form of derivatives of the principal Jacobian 
J. Dealing with the generalized stream function £ = {£ , £1, £21 £3}) the following 
identities are useful 

9 ^o, ftgjg^O (2.4) 

d 2 J 1 f dJ dJ dJ dJ \ , . 

J ( "n7 "n7 "o7 ) (2-5) 



See details of working with Jacobians and the generalized stream functions in [H 
Rylov,2004]. 

Example. Application of the stream function for integration of equations, describ- 
ing the 2D stationary flow of incompressible fluid. 
Dynamic equations have the form 

U X +Vy = 0, 8y (UU X + VUy) = 8 X (UV X + VVy) (2.6) 

where u and v are velocity components along x-axis and y-axis respectively. 
Introducing the stream function if) by means of relations 

u = -il) v , v = ip x (2.7) 

we satisfy the first equation (12.61) identically, and we obtain for the second equation 
(12.61) the relations 



*Py (iPxyy + ^ xxx) = ^x {^xxy + 4 >yyy) 



which can be rewritten in the form 

d(x,y) 



0, u = ^ xx +i/; n (2.8) 



where to is the vorticity of the fluid flow. The general solution of equation (I2.8P has 
the form 

w = v„ + K = n W (2.9) 

where Q is an arbitrary function of i/j. 

For the irrotational flow the vorticity Q (ip) = 0, and we obtain instead (12. 9p 

V>xx + Vw = o (2.10) 

One obtains the unique solution of (I2.10p inside of a closed region of 2D space 
provided, that the value "0|e of the stream function ip is given on the boundary S of 
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this region. The differential structure of equations (12.91) and (12.101) is similar. One 
should expect, that giving the value of the stream function ip on the boundary E, 
one obtains the unique solution of the equation ( 12.101) . But it is not so, because the 
indefinite function Q (ijj) is not given, and it cannot be determined from the boundary 
condition, because the nature of the function Q (-0) is another, than the nature of 
the boundary conditions. First, if the flow contains closed stream lines, which do 
not cross the boundary, one cannot determine the values of Q on these stream 
lines from the boundary conditions. But for determination of the unique solution 
the values of Q on the closed stream lines must be given. Second, the boundary 
conditions are given arbitrarily. The function Q cannot be given arbitrarily. For 
those stream lines, which cross the boundary more than once, the values of fl on 
the different segments of the boundary are to be agreed. Thus, the nonuniqueness 
of the solution, connected with the indefinite function Q has another nature, than 
the nonuniqueness, connected with the insufficiency of the boundary conditions. 

3 Derivation of hydrodynamic equations from the 
variational principle 

We use the variational principle for the derivation of the hydrodynamic equations 
(11.11) . (11.21) . (jl.4p . The action functional has the form 

A[U,P] = J {f p -P E ^-P k (j k -P^^^)} d ^ ^ 

where pf., k = 0, 1, 2, 3 are the Lagrange multipliers, which introduce the designations 
for the 4-flux 

f) j 

3 k = Po(Z)j^, A; = 0,1, 2, 3 (3.2) 

Note, the expression for the 4-flux (13.21) satisfies the first equation (12. ip identically, 
because the expression (I3.2p may be reduced to the form of the second relation (12. ip 
by means of a change of variables £ 

6o = 6o, 6i = J Poi€)d£i, 62 = 62, 63 = 63 

Then 

d (60,61, 62,63) d (60,61,62,63) ~ r 

Po (£) a 1 i 2 3T = a ( n 1 2 r\"> 61 = Po {£) d£i 
o [x u , x 1 , x z , x 6 ) a (x u , x , x z , x 6 ) J 

Besides according to the first identity (12.41) the relation (13.21) satisfies the Lin 
constraint (ll.4p . 

Variation of the action (13.11) with respect to pk = {po,p} gives relations (13.21) . 
Another dynamic equations have the form 

s »-- v° = -h-T P {pE(p)) = -^-T P (pE) (3 ' 3) 
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-d. (po «)ft /l ) +P*^T (0 #- = °> « = 0,1,2,3 (3.5) 
Using the third relation (12.41) . we obtain 

#J ^"'^^Mfl^+P.llfl^-O (3.6) 



- J 



Now using (12.51) . we obtain 

dJ dJ dJ dJ \ fdp (£) 



-i 



°' fe (3.7) 



Using the first relation (|2.4p . we obtain 

■ r, (£S-£S ft,(ete = ° 1 ( = 0,1,2,3 (3-S) 

There are two ways of dealing with this equation: 

1. Elimination of GSF £, which leads to the Euler equations. 

2. Integration, which leads to appearance of arbitrary functions. 

The first way: elimination of GSF 

Convoluting (13. 8p with and using dynamic equations (13.21) . we obtain 

j k d iPk -j k d kPi = 0, 1 = 0,1,2,3 (3.9) 

Substituting p k and j k from relations H3. 31) and (13 .411 . we obtain the Euler dy- 
namic equations (11.11) 

f) 1 3W 

d v a + (vV)v a = -d a —(pE) = —d a p, a = 1,2, 3, p = p 2 — (3.10) 

op p op 

The continuity equation (II. 2p is a corollary of equations (13.21) and identity (12.1 1) . 
Finally the Lin constraints (11.41) are corollaries of the first identity (12. 4p and dynamic 
equations (13. 2p . 

The second way: integration of the equation for p k 

Let us consider the equations (13.81) as linear differential equations for p k . The 
general solution of (13.81) has the form 

Pk = (dk<p + g a (£) d k U , k = 0, 1, 2, 3 (3.11) 

where g a (£) , a — 1,2,3 are arbitrary functions of £, (p — g° (£ ) is a new vari- 
able instead of fictitious variable £ . Let us differentiate (13.111) and substitute the 
obtained expressions 



d sPk = (d s d k <f + g a (£) d s d k U + ^Md^d^ (3.12) 



in (13.81) . Using the first identity (12.41) . we see, that the relations (13.121) satisfy the 
equations (13. 8p identically. 

We may substitute (13.111) in the action (13.11) . or introduce ( 13. lip by means of the 
Lagrange multipliers, (the result is the same). We obtain the new action functional 



L- pE (p) - j k [d W + g a {i) d k U J d\ (3.13) 



which contains arbitrary integration functions g (£). Here 

f = P, j = pv = {j 1 ,j 2 ,j 3 } (3.14) 
The integration functions g (£) are considered as a fixed functions of The term 



PkT^— = (<9 fc ^ + g (£) <9 fe £J 



0,fe 



d €o.k d(x Q ,x l ,x 2 ,x 3 ) 



is omitted, because it does not contribute to dynamic equations. 
Variation of (13 . 1 3[) with respect to (p, p and j M gives respectively 



5<P ■ d k j k = 



Sp : d p + g? (€) d ^ + ^ + ^-(pE(p))=0 

W ■ v^ J j = d^ + g a (0 d,i a 
Variation of (13 . 1 3[) with respect to £ a gives 

5L ■ (o m a + (vv) a = o, 

where 



dg a (£) 00" (£) 



and v is determined by the relation ( 13. 18ft 
If det I If2 a/3 | I ^ 0, then the Lin constraints 



(doL + (W) U = 

follows from (13.191) 

However, the matrix Q al3 is antisymmetric and 



det 



a/3 I 






n 12 


Q 13 


-n 12 





n 23 


-n 13 


-n 23 






(3.15) 



(3.16) 
(3.17) 
(3.18) 

(3.19) 
(3.20) 



(3.21) 



(3.22) 
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Then it follows from (13.191) 



doL + (W) C = r^Sa^ (0 a = 1, 2, 3 (3.23) 

where a; = a; (t, £) is an arbitrary quantity, and p (£) is the weight function from 
(E2). 

The obtained equation (13.231) contains the initial dynamic equation (11.41) as a 
special case. For irrotational flow, when VL^ 1 (£) = 0, the equation (13.231) turns 
to (jl.4p . In the action functional (13. ip the initial relation (11.41) is used as a side 
constraint. It is a reason, why the equation (I3.23P is not obtained from the action 
functional (13.11) . 

Note, that eliminating the variables <p and £ from dynamic equations (13 . 1 71) - 
(I3.19p . we obtain the Euler dynamic equations (11.11) . 

The vorticity w = V x v and v x uj are obtained from (I3.18P in the form 

u;o = V x v = x V£ Q (3.24) 

v x cv = ^V^(vV)C (3.25) 

Let us form a difference between the time derivative of (13.181) and the gradient of 
(I3.17p . Eliminating Q afl (£) do$, a from the obtained equation by means of equations 
(I3.19p . one obtains 

<%v + Vy + Vp - ^V^vV)^ = o (3.26) 

Using (I3.25P and (I3.24p . the expression (13.261) reduces to 

v 2 1 OE 

9 v + V— + -V(p 2 — )-vx (Vx v) = (3.27) 
2 p op 

In virtue of the identity 

v 2 

v x (V x v) = Vy - (vV)v (3.28) 

the equation (I3.27P is equivalent to (II. ip . 

Note, that the Euler equations (II. ip are obtained at any form of the arbitrary 
function uj (t, £) in the equations (13.231) . because the equations (13.231) are used in the 
form (13.191) . where the form of uo (t, £) is unessential. Solution of the Euler system 
(11.11) . (11.21) in the form p = p (t, x), v = v (t, x) does not depend on the form of the 
indefinite function uj (t, £). 

If oj (t, £) ^ 0, the dynamic equations (13.231) describe a violation of the Lin 
constraints (11.41) . One obtains another labeling of the stream lines, than that one, 
which is described by the Lin constraints (11.41) . If the flow is irrotational, and ft = 0, 
the labeling does not depend on u (t, £). 
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Let us consider two different labeling £ and | of the the same fluid flow described 
by the variables p = p(t, x), v = v(£, x). The initial conditions are supposed to 
have the form 

<p(0,x) = <p in (x) = 0, p(0,x)=p in (x), (3.29) 
L(0,x) = (L) in (x)=a; a , a = 1,2, 3 (3.30) 

Then according to (13TT8D . (13T20D 

# Q (£ in (x)) = ^ (x) = < (x) , a = 1, 2, 3 (3.31) 

According to (13.231) the dynamic equations for labeling £ and £ have the form 

<9oC + (WK Q = « = 1,2, 3 (3.32) 

dol a + (W) £ Q = ^rs api n^ (|) a = 1, 2, 3 (3.33) 

Po I? J 

where the velocity v = v (t, x) is the same in both equations and the function Vt" 1 (x) 
is defined by the relation 

(3.34) 

If the velocity v is defined by relations (I3.18p . it satisfies the Euler equations and 
associates with the generalized stream function £ (t, x), whose evolution is described 
by the equations (I3.23f) 

In general, the evolution of the quantities £ and £ is different, although the 
coincide at t — 0. Let 

r/ = f-£ (3.35) 

It follows from (13.321) and (13.331) that mismatch r\ between \ and £ is determined 
by the relation 

d oVa + v (t, x) V Va + - ^ gq^Qfr (| (t, x)) = 0, r? Q (0, x) =(8,36) 

a = 1,2,3 

The system of ordinary differential equations, associated with the equation (13.361) . 
has the form 

dt dx 1 = dx 2 = dx 3 = Po (I (*,x)) cfr^ g7 

1 ^(i,x) ^ 2 (t,x) t;3(t,x) w(t,£(t,x))e^fi^(£(t,x))* ' 

a = 1,2,3 

Solution of the system of ordinary equations at the initial conditions rj (0, x) = 
has the form 

t 

r ^(^M^ ^ (| (t> x)) a, ^ = 1,2,3 (3.38) 
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Then 



t 

t» (t, x) = ^ (t, x) - / ^^^ x ))^ ^7 (| (t) x) ) dt> M =l,2,3 (3.39) 

Thus, although solution p, v of the Cauchy problem for the Euler system of 
hydrodynamic equation (11.11) . (jl.2p is unique, the solution p, v, £ for the Cauchy 
problem of the complete system of hydrodynamic equations (jl.ip . (11.21) . (I3.23P is 
not unique. The reason of this nonuniqueness is consideration of interfusion. This 
consideration is formal. One cannot understand mechanism of the interfusion influ- 
ence from this consideration. Nevertheless this influence takes place, and it should 
be investigated more closely. 



4 Two-dimensional vortical flow of ideal barotropic 
fluid in the three-dimensional space. 



(4.1) 



It seems, that in the two-dimensional flow instead of determinant (I3.22p we have the 
determinant 

o n 12 
-n 12 o 

which does not vanish, in general. Then the problem of nonuniqueness of thelabeling 
is removed and the solution of the Cauchy problem for the complete hydrodynamic 
system becomes to be unique. 

In reality, we may control the solution only via initial conditions. We may give 
the two-dimensional initial conditions, i.e. 



<9 3 v in (x) = 0, <9 3 Pin ( x ) = 0, 



vl (x) 



In this case 



K (0 = o, 



The determinant 



det\\n af3 \ 



Q 12 

-n 12 o o 





= 



and the relations ( 13.23P take the form 

^!+v(t,x)V^ = 0, 

9 e 3 + v(t,x) v^ 3 



<9 £ 2 + v(t,x) V£ 2 = 



(4.2) 
fi£(*)=0 (4.3) 



(4.4) 

(4.5) 
(4.6) 



' 3 Pq(£i,£ 2 ) 

One cannot control indefinite quantity u (t, £), which may depend on x 3 . The equa- 
tion (14. 6p generates nonunique solution of the Cauchy problem of vortical flow for 
the complete hydrodynamic system. The flow with the two-dimensional initial con- 
ditions turns into three-dimensional vortical flow. 
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5 Concluding remarks 



Solution of the Cauchy problem for the vortical flow of inviscid barotropic fluid is not 
unique, if we solve seven dynamic equations of the complete hydrodynamic system, 
which includes description of the shape of stream lines and has seven dependent 
variables p, v, Nonuniqueness is connected with the fact, the initial conditions 
for variables p, v, £ do not control the intermixing effect. 

Solution of the Cauchy problem for the vortical flow of inviscid barotropic fluid 
is unique, if we solve only four dynamic equations of the Euler system and ignore 
shape of stream lines. In this case dynamic equations are written for four dependent 
variables p, v. 

Solution of the Cauchy problem for the irrotational flow of inviscid barotropic 
fluid is unique for both seven equations of the complete hydrodynamic system and 
for four equations of the Euler system. 

The intermixing effect, generating nonunique solutions, associates with the tur- 
bulence phenomenon at the following points: (1) both effects are not controlled by 
the initial data for variables p, v, £ of the hydrodynamic equations in the Euler 
representation, (2) both effects take place at the vortical flows, and they are absent 
at the irrotational flows, (3) both effects are strong at vanishing viscosity. We admit, 
that the interfusion may be connected with the turbulence phenomena, although we 
do not yet insist on this statement. 

If, indeed, the turbulent phenomena are connected with the interfusion and with 
the shape of stream lines, it becomes clear, why numerous investigations of hydrody- 
namic equations, describing only density and velocity, but not the shape of stream 
lines, had not led to a progress. The researchers looked for turbulence in that region, 
where it is not placed. 

References 

[1] A. Clebsch, Uber eine allgemaine Transformation der hydrodynamischen Gle- 
ichungen, J. reine angew. Math. 54 , 293-312 (1857). 

[2] A. Clebsch, Ueber die Integration der hydrodynamischen Gleichungen, J. reine 
angew. Math. 56 , 1-10, (1859). 

[3] C.C. Lin, Hydrodynamics of Helium II. Proc. Int. Sch Phys. Course XXI, pp. 
93-146, New York, Academic, 1963. 

[4] Yu. A. Rylov, Hydrodynamic equations for incompressible inviscid fluid 

in terms of generalized stream function. Int. J. Math. & Mat. Sci. 

vol. 2004, No. 11, 21 February 2004, pp. 541-570. (Available at 
http: / / arXiv.org/abs/physics/0303065 ) . 

[5] R. Salmon, Hamilton fluid mechanics, Ann. Rev. Fluid. Mech. 20, 225-256, 1988. 



12 



CAPTIONS 

Figure 1. Dashed lines show real trajectories of particles. The solid lines show 
trajectories of the mean particle motion. 
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Fig. 1 . Two noninteracting flows. Dushed 
lines show real motion of particles. Solid 
lines show mean motion of particles (lines 
of constant labelling) 



